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We develop a semiclassical method for the determination of the nonlinear dynamics of dissipative 
quantum optical systems in the limit of large number of photons N, based on the 1/iV-expansion and 
the quantum-classical correspondence. The method has been used to tackle two problems: to study 
the dynamics of nonclassical state generation in higher-order anharmonic dissipative oscillators and 
to establish the difference between the quantum and classical dynamics of the second-harmonic 
0^ ' generation in a self-pulsing regime. In addressing the first problem, we have obtained an explicit 

time dependence of the squeezing and the Fano factor for an arbitrary degree of anharmonism in the 
short-time approximation. For the second problem, we have established analytically a characteristic 
time scale when the quantum dynamics differs insignificantly from the classical one. 
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I. INTRODUCTION 

> ' 

The situation when nonlinear interactions involve a large number of photons, N, is quite typical of many of the 
problems of quantum and nonlinear optics Jl]-||]. Hiedmann et al. suggested [Q to use the method of the 1/iV- 
expansion |5) to describe the nonlinear dynamics of the mean values and second-order cumulants of a quantum system 
in the limit N 3> 1. Following the general scheme of that method an exact or approximate solution can be found 
' first in terms of the coherent state representation in the classical limit N — * oo and then adjusted by adding the 
. quantum corrections. The method proves to be particularly convenient when the generation dynamics of nonclassical 
states needs to be determined Q. We have recently developed the method further to study the enhanced squeezing 
at the transition to quantum chaos 

Papers are concerned with consideration of the problems of nondissipative quantum systems only. In this 

paper we extend the method to dissipative quantum systems. For quantum systems without dissipation, the lowest- 
order of the l/A^-expansion is equivalent to the linearization in terms of the classical solution whereas in 
dissipative systems, as will be demonstrated herein, the solution of motion equations for variations near the classical 
' trajectory cannot provide complete information on the dynamics of quantum fluctuations even in the lowest order of 
^ . l/N. We will show that the influence of reservoir on the dynamics of the expectation values and dispersions, which 
K^J ' is different from the energy dissipation, always exists; it has the quantum nature and can not be neglected even in 
the semiclassical limit. Specific manifestations of the effect will however depend on the type of the attractor in the 
underlying classical dynamic system. For systems with a simple attractor in the classical limit, the "quantum diffusion" 
associated with the reservoir quantum fluctuations will not introduce any new physical effects in the dynamics of the 
main system, at least in the short-time limit. In contrast to that, for a stable limit cycle, such a diffusion appears to 
be the main mechanism responsible for the difference between the classical and quantum dynamics at N 3> 1. 

Along with the presentation of a general formalism, two typical examples of quantum optical systems with a simple 
attractor and a stable limit cycle in the classical limit N — > oo will be considered: the dissipative higher-order 
anharmonic oscillator and the self-pulsing regime of intracavity second-harmonic generation (SHG). We will show 
how the method of 1/A^-expansion can be used to investigate the dynamics of nonclassical state generation and to 
determine a time scale for a correct classical description of the dissipative quantum dynamics. 

A quantum anharmonic oscillator with the Kerr-type nonlinearity is one of the simplest and most popular models 
employed in the description of quantum statistical properties of the light interacting with a nonlinear medium [QU ■ 
The Kerr oscillator model with a third-order nonlinearity yields an exact solution in both nondissipative and 
dissipative limits ||. However, due to the complexity of the solution in the dissipative case, numerical methods or 
special approximate analytical methods have to be used to determine statistical properties of the radiation in a most 
relevant experimental case involving a large number of photons. Moreover, there are no exact solutions available for 
the model of the anharmonic oscillator with a higher-order nonlinearity. 

In this paper, a simple and explicit time dependencies of the degree of squeezing and the Fano factor are obtained 
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analytically in the model of anharmonic oscillator of an arbitrary order for the most interesting experimental situation 
featuring higher intensities (N 3> 1) and short-time interactions. As another example of application of the l/N- 
expansion, we consider the self-pulsing in SHG O]. Such an oscillatory regime, corresponding to the limit cycle, was 
observed experimentally in fll2| . There are quite a few papers dealing with the development of approximate analytical 
and numerical methods with the purpose of describing different dynamic regimes in SHG in terms of quantum 
mechanics [p~3|— p~T| . In particular, Savage O] calculated the Q distribution function in the Gaussian approximation 
about the classical limit cycle. He demonstrated numerically that in the classical limit, the initial rapid collapse of the 
Q distribution in the neighbourhood of the limit cycle was followed by the diffusion around the limit cycle. However, 
the author did not offer any analytical solution of the problem or explanation of the physics of the effect observed. 

In this paper we show that the diffusion around the classical limit cycle can be obtained as a solution of the motion 
equations for low-order cumulants by using the 1 /TV-expansion technique. This enables us to find the time scale 
t <C f with t* ~ 2A r 7~ 1 (7 is a damping constant), for a correct classical description of self-oscillations in SHG. 
The resultant estimate is consistent with that obtained for t* numerically in Finally, we interpret the quantum 
diffusion around the limit cycle as having been caused by the effect of the reservoir vacuum on the SHG dynamics. 

The structure of the paper is as follows. Section || describes a general formalism of the 1/iV-expansion applicable 
to an arbitrary single-mode quantum dissipative system and presents the solution of the motion equations for mean 
values and second-order cumulants obtained in the first order of 1/N. Sections III and [fy| deal with the nonclassical 
state generation dynamics in higher-order anharmonic oscillators and the quantum-classical correspondence for the 
self-pulsing regime in SHG, respectively. The final section contains a summary and concluding remarks. 



II. 1 /TV-EXPANSION AND QUANTUM-CLASSICAL CORRESPONDENCE 

First of all, we need to generalize the approach of [Q for the case of systems with dissipation. As an illustrative 
example we consider a quantum anharmonic oscillator with the Hamiltonian in the interaction picture 

H = Atfb+j^ [ (tfb) l+1 , [b,tf]=l, (1) 

where the operators b and describe a single mode of a quantum field and the constant A; is proportional to a 
(2/ + l)-order nonlinear susceptibility of a nonlinear medium (I is an integer), A is the light frequency detuning from 
the characteristic frequency of quantum transition, and h = 1. Everywhere in this paper we use the normal ordering 
of operators. The oscillator interacts with an infinite linear reservoir of finite temperature. The Hamiltonians of the 
reservoir and interaction of the oscillator with reservoir are defined as follows 

H r = ]T ^ {d)dj + 1/2), H int = M^ f + H.c.) , (2) 

3 3 

where the Bose operator dj ([dj, d\.] — Sjk) describes an infinite reservoir with characteristic frequencies ipj, and Kj 
are the coupling constants between reservoir modes and the oscillator. Introduce new scaled operators a = b/N 1 / 2 , 
cj = dj/N 1 / 2 and the Hermitian conjugates satisfying the commutation relations 

[a,J] = l/N, [ Cj ,4]=6 jk /N. (3) 

In the classical limit N — > 00, we have commuting classical c- numbers instead of operators. Now the full Hamiltonian 
H = H + H r + Hint m ay be rewritten as H = NTC, where TL has the same form as (|l|) and (||) but for the following 
replacements 

b^a, tf^a f , dj^cj, d] -*■ c), and A; -> gi(N) = XN l . (4) 

It can be shown that the photon- number dependent constant gi(N) provides a correct time scale of oscillations for the 
nonlinear oscillator (|l|) in the classical limit (for the case of Kerr nonlinearity with 1 = 1, see, e. g. [Q). Note that Ti. 
can have an explicit time dependence in the general case . Within a standard Heisenberg-Langevin approach, the 
equation of motion has the form ( chap. 7) 

a = -i(A-i^)a + V + L(t), (5) 
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where V = dHo/da^, 7 = 2it\k(uj)\ 2 p(lu) is the damping constant, p{u>) being the density function of reservoir 
oscillators, which spectrum is considered to be flat. The Langevin force operator L(t) is in a standard relation to the 
operators {cj} of the reservoir |Q. The properties of L(t) |lj in our notations (||) may be rewritten as 

(£(«)>« = (tf(t)) R = 0, (i t a) fi +(at£} R = 7 -^, (La) fl + (aL)* = 0, (6) 

where the averaging is performed over the reservoir variables and (rid) is a single-mode mean number of reservoir 
quanta (phonons), that is related to temperature T as (rid) — [exp (w) — l] j where k is the Boltzmann constant 
and u) is the characteristic phonon frequency. From the Heisenberg-Langevin equations for a, a 2 and the Hermitian 
conjugated equations, by using Eqs. (||) and (||), we obtain: 

i±(a) = { V)-il(a), 

l i {{Sa)2) = 2{v6a) + {w) ~ i ^ 5a ?^ ( 7 ) 

i^(Sa*Sa) = -(V*Sa) + (5a*V) - ij(Sa*6a) + ij^ c 



dt 1 ' ' ' JV ' 

where W = (l/N)dV/da\ z = (a), (((5a) 2 ) = (a 2 ) — z 2 , (5a*<5a) = (aJa) — \z\ 2 , and the averaging is performed over 
both the reservoir variables and the coherent state |a) = exp(Naa' f — 7Va*a)|0) corresponding to the mean photon 
number ~ N. In deriving Eq. (^), we neglect the insignificant additional detuning introduced to A by the interaction 
with the reservoir 0. In the absence of no damping 7 = 0, our equations for the mean values and the second-order 
cumulants ^ are reduced to the corresponding equations in j|,f[| . 

Set of equations (Q) is not closed and is basically equivalent to the infinite dynamical hierarchy system for the 
cumulants of a different order. To truncate it up to the second-order cumulants, we make the substitution a — > z + Sa, 
where at least initially the mean z ~ 1 and the quantum correction \Sa(t — 0)| ~ N~ x / 2 -C 1. Using the Taylor 
expansion of the functions V and W and after some algebra analogous to that used in , we get from ^ in the first 
order of 1 /N the following self-consistent system of equations for the mean value and the second order cumulants (for 
details see [@) 

iz = -ilz + (V) z + ±Q(z, z\ C, C, B), (8a) 



and the corresponding equation for C*(t) that could be obtained from equation (Sb) by way of complex conjugation 



The quantum correction to the classical motion Q in Eq. (8a) has the following form 



K£). c+ K£).H^).H''- 



In Eqs. (g) and (g) the subscript z means that the values of V and its derivatives are calculated for the mean value 
z and we have introduced 

B = N(5a*5a) + 1/2, C = N((Sa) 2 ). (10) 
The initial conditions for system (H) are 

5(0) = 1/2, C(0)=0, (11) 



and an arbitrary z(0) = zq which is of order unity. The equilibrium value of cumulant B in equation ( |8c[ ) is determined 
by the mean number of reservoir's quanta and its zero-point energy as 
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£?<°> = (n d ) + 1/2. 



(12) 



Note that the zero-point energy of the reservoir appears in the equations of motion for the cumulants though it was 
not presented in the Heisenberg equations of motion and even may be dropped from the Hamiltonian redefining a zero 
of energy. Such "reappearance" of a zero-point field energy is quite common in other problems of quantum theory 
where a vacuum is responsible for the physical effects [M . 

The motion equations for second-order cumulants B and C [Eqs. (8b), (pq)] are linear inhomogeneous equations. 



Their solution consists of two parts: a general solution of the homogeneous set of equations (i.e. without term +ijB^ 
m Eq. @) that we denote as (B(t), C(*)J, and the particular solution of the inhomogeneous equations 

(B(t), C(t)) = (B(t)Mt)) + (lB(°h, 0) . (13) 

To find (B(t), C(t)) we use the perturbation theory for N ^> 1 and as a first step neglect the quantum correction 
Q/N in Eq. (|8a|). It is easy to see that the homogeneous equations of motion for cumulants ( jSq) and (jSc)) can be 
obtained from the classical equation (i.e. from (|8a| ) with Q/N — > 0) by linearization around z (substitution z — > z + Sz, 
\Sz\ <C \z\), if one writes the dynamical equations for the variables (5z) 2 and |(5zP. The only difference between the 
linearization of classical motion equations and equations for quantum cumulants ( |8b[ ) , ( |8c|) lies in the impossibility to 
get the initial conditions ([ll]) for C and B from only initial conditions for linearized classical equations of motion (see 
also the discussion of this problem in 0). Hence, we first need to know the classical solution z c \(t), find differentials 
dz c \ and dz* x , and then use the substitution (B(t),C(t)) -> (|dz| 2 , (dz) 2 ). 

Thus, it has become apparent that assuming the actual field deviates little from the coherent state and treating 
the small deviation as as a first-order correction would not be equivalent to direct linearization around a classical 
trajectory. Even in the limit N — > oo, we will always deal with the influence of reservoir on the dynamics of the 
quantum system via the second-order cumulant B, which has the form of quantum diffusion 

B(t) = B(t) + «n„) + 1/2)7*, (14) 

where B is obtained from linearization around a large mean field. In particular, as follows from Eq. (|l4|), the quantum 
diffusion also exists for the case of a quiet reservoir (no) = 0. 

We now discuss the range of validity of the 1 /iV-expansion and the role of quantum diffusion in different classical 
dynamical regimes. The criterion of validity of the l/7V-expansion may be represented in two forms. First, the 
1/iV-expansion works well, provided the difference between the classical and quantum solutions is small 



z(t) - z c \(t) 



z c \(t) 



1 



/* Q(t')dt' 



N \z(t)\ 



« 1, (15) 



where z c \(t) is the solution of Eq. (pa) for N — > oo. To write the second form of the criterion of validity of the 

1 /TV-expansion we introduce following ||,[7j the "convergence radius" R = {[Re(5a)} 2 + [lm(Sa)} 2 } 1/2 . Then, the 
expansion is correct within a time interval when 

R(t) B^ 2 (t) 



\z(t)\ NV*\z(t)\ 



< 1. (16) 



As_a rule, the both conditions, Eqs. ( |lq ) and ([Tq), determine the same time interval for validity of the 1 /TV-expansion 
(For a physically interesting exception, the problem of SHG, see Sec. [Tv| ). 
For dissipative systems with a simple attractor, the classical field intensity |z c i(*)| 2 , as well as cumulants B(t), C(t) 
and quantum correction Q(t ) are proportional to the factor exp(— 7*) and therefore, as follows from Eqs. (JlJ) and 



( |16| ) with account of Eq. (14), the 1/iV-expansion is well defined only in the time interval of order of several relaxation 
times: t* ~ 7 _1 p9[ . Moreover, during this time interval the influence of quantum diffusion on the system dynamics 
is small. 

A quite different behavior is characteristic for the stable limit cycle. Here, a variation near classical trajectory 
collapses to zero (5a — > 0) and therefore B(t) ~ \Sa\ 2 — > 0, C(t) ~ (5a) 2 — > 0. However |z c i(*)| — 1 for the limit cycle 
and, as a result, the time interval of validity of the l/AT-expansion is fairly large ** ~ N^ 1 . What is important, the 
diffusion is the major physical mechanism responsible for the difference between the classical and quantum dynamics 
for a stable limit cycle. In the two following sections we consider two typical examples of dissipative optical systems 
with a simple attractor and limit cycle. 
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III. NONCLASSICAL STATES GENERATION IN HIGHER-ORDER ANHARMONIC OSCILLATORS 



We start by defining the squeezing and the Fano factor. Define the general field quadrature as Xg — aex-p(-iO) + 
cv exp(i9), where 8 is the local oscillator phase. A state is said to be squeezed if there exists some value for which 
the variance of Xg is smaller than the variance for a coherent state or the vacuum jl]]9| . Minimizing the variance of 
Xg over 0, we get the condition of so-called principal squeezing ]l|,P,|To| in the form 

S^l + 2N((\Sa\ 2 ) - \{(5a) 2 )\) = 2(B \C\) < 1. (17) 

The determination of the principal squeezing S is very useful because it gives the maximal squeezing measurable by 
the homodyne detection [|||. 

Another important characteristic of nonclassical properties of the light is the Fano factor F — ((n 2 ) — (n) 2 )/(n), 
that determines the deviation of probability distribution from the Poissonian . Substituting expression a — ► z + 8a 
into the expressions for (n) = N(a'a) and (n 2 ) = N 2 (ctfaa*a) = N 2 (a? 2 a 2 ) + (n) and after the Talor expansions in 
the first order of 1/N, we have 

F = 2B + ( — C + c.c. ) . (18) 



We see, that in order to determine the time dependence of the principal squeezing S in ( |17| ) and the Fano factor 
( |l8| ) for nonlinear oscillators, we need to find the time dependence of z, C, and B from (|8|) for the Hamiltonian (|l). 
Following the general procedure described in previous section, we first neglect quantum correction Q/N in Eq. (|8a). 
In this case, equation (|8a|) has an exact solution in the form 

z(t) = ^ cxp [HA - 7 /2)*] exp [-i gi \z f l m(t)] , /i/(*) = [1 - exp(- 7 Zt)] / 7 L (19) 

We find the differentials dz and dz* of classical solution (p9|), and using the substitution \dz\ 2 + B — > B and (dz) 2 — > C, 
we get 

C(t) = -lz 2 \z \ 2 ^ gm (t) (l\z \ 2l g ltM (t) + i) exp [(- 7 - i2A)t - i2\z \ 21 gm (t)] , 

B(t) - exp(- 7 *) [1/2 + l 2 \z fg 2 tf (*)] + «n„) + 1/2) yt, (20) 



where we have taken into account the initial conditions for B and C, Eq. (11). By substituting formulas (^0|) into 
Eq. (jl7|), we obtain in the limits r = gi(N)t -C 1 and 7 * -c 1 a very simple dependence of S on time as 

S(t) = 1 - [lx 21 - (l/gi)(n d )] 2t < 1, (21) 

where for the sake of simplicity we have assumed that the initial value zo is real, xo = Rezo, and we have taken 
into account only terms that are linear in r and 7*. The short-time approximation r -C 1 as well as the limit of 
a large photon number N 3> 1 are quite realistic for a nonlinear medium modelled by the anharmonic oscillators 
(for numerical estimates, see Q, chap. 10, and jl0|]). It should be noted that our formula ( |2l| ) coincides with the 
corresponding formula for S(t) in jlO) for the Kerr nonlinearity (I = 1) without loss (7 = 0). In the case of no loss 
(7 = 0), our formula ( |2l| ) shows that the rate of squeezing is determined by the factor 21xq 1 XiN 1 = 21V < ' 21+1 ^ . Since 
A; is proportional to the (21 + l)-order nonlinear susceptibility, the factor T>^ 21+1 ^ has a physical meaning of nonlinear 
polarization. Therefore, the stronger is nonlinear the polarization induced by light in the medium, the more effective 
squeezing of light is possible. For a finite dissipation 7 / 0, the squeezing is determined by an interplay between the 
polarization of nonlinear medium modelled by the anharmonic oscillator and the thermal fluctuations of the reservoir. 
As follows from (pl|), there exists the a critical number of phonons (rtd)^' = (l/^V^ 21 ^ such that for (rid) > ("■d)^ cr ' ) 
the squeezing is no longer possible. 

In the same approximation, we obtain from (h8[) the following time dependence of the Fano factor 



F(t) = 1 + 2(n d )7*. (22) 

Thus, the statistic is super- Poissonian for any 7^0 and is independent of the degree of nonlinearity I. This is 
in a good agreement with the earliest result of Q for the case of a dissipative Kerr oscillator (I = 1), where the 
impossibility of sub-Poissonian statistics and antibunching were found from the exact solution. 

Turn now to the discussion of the ranges for validity of our approach. It is easy to see that in terms of our approach 
the time dependence of the number of quanta for / = 1 is 
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(n)(t) + 1/2 = N\z\ 2 + B^N\z \'(l- 1 t) + (71^1, 7 t « 1, 9l t < 1, (23) 



where we have used expressions (|20| ) for cumulants B and C. It is instructive to compare fl23| ) with the exact solution 
for (n) (t) for the Kerr nonlinearity M 

(n)(t) - (n ) exp(-Tt) + [1 - exp(- 7 i)](n d ). (24) 

Eq. (^3|) and Eq. (^4|) both describe the evolution of an initially coherent state to a final chaotic state being 
characteristic for the reservoir. It is evident that formulas ( |24] ) and (^) coincide, when yt <C 1 and (no) — iV ^ 1. A 
more accurate analysis of the condition for validity of the l/7V-expansion should include a comparison of the solution of 
quantum motion equation ( |8a|), w hich takes into account the quantum correction Q/N given by (|9|), with the solution 
of classical motion equation (|19|). It may be shown after some algebra, that if yt <C 1 and r <C 1, the influence of 
the quantum correction Q/N on the dynamics of the mean value z is of the order 1/N and, therefore, our cumulant 
expansion is well-defined for N ^> 1. The same conclusion could be obtained considering another criterion of validity 

IV. QUANTUM-CLASSICAL CORRESPONDENCE IN SELF-PULSING REGIME OF 

SECOND-HARMONIC GENERATION 

We now consider another example of a quantum optical system, namely intracavity SHG. The Hamiltonian describ- 
ing two interacting quantum modes in the interaction picture has the form |Tl| , |l4|| 

H = J2 A 3 -&& + iEN^{b\ h) + |(&f& 2 b\b\), (25) 
i=i 

where the boson operators bj (j = 1,2) describe fundamental and second-harmonic modes, respectively, Aj is the 
cavity detuning of mode j, EN 1 / 2 is the classical field driving first mode (E is of order of unity), x is a second-order 
nonlinear susceptibility. The linear reservoir and its interaction with a second-order nonlinear medium are described 
by the Hamiltonians (||). Now we can rewrite full Hamiltonian of the problem in the form H = NH, where TL has the 
same form as ( p5| ) and (Q) with account of the replacements analogous to (^) and definition of new coupling constant 



9 = XVN, (26) 

which is of order unity. Formally, the procedure of the 1/iV-expansion developed in sec. || can not be applied to 
the problem of SHG, however its straightforward generalization to the case of two interacting modes gives in the first 
order of 1 /N the following self-consistent set of equations 

z 1 = -'^-z 1 +E + gz* 1 z 2 + -^gB 12 , (27a) 

72 9 2 1 9 n / _, s 

Z2 = --z 2 --z 1 ---C 1 , (27b) 

Bx - - 7l (B a - BC°)) + g B* l2 zi + gBml + C{z 2 + d**, (27c) 

B 2 = - 12 {B 2 - fl(°)) - gB\ 2 z x - gB 12 zl, (27d) 

Ci = -yid + 2g(C 12 z* 1 + B x z 2 ), (27e) 

C 2 = - l2 C 2 -2gC 12Zl , (27f) 

Cjs = -0.5(7i + 72 )Ci2 + gB 12 z 2 - Qzi + C 2 z{, (27g) 
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B12 = -0.5(71 + j 2 )B 12 + gC 12 z* 2 +gz 1 {B 2 - Bi), (27h) 



where Zj = (a 3 ) = N 1 / 2 ^}, Bj = N{8a* j 8a j ) + 0.5, C s = iV((5oy) ) (j = 1,2), S i2 = iV^aJfas), C12 = N(S ai 6a 2 ) 
and is denned in Eq. (|l§). The initial conditions for system @ are Bj(0) = 1/2, (7,-(0) = C 12 (0) = B 12 (0) = 0, 
^2(0) = 0, and zi(0) = zq, where z$ is of order unity. In this work, we limit ourselves only by the values of the field 
strength z a corresponding to self-oscillations Q and A x = A 2 = 0. 

It is easy to see that in the limit N — > 00 and for g — const ~ 1, we get from Eqs. ( 27a| ) and (27b) correct classical 



motion equations for the scaled field amplitudes. The solution of motion equations (27c)-(27h) for the second-order 
cumulants has the form 

X(i) = X(t) + ( 7 £^i,TB(°^,0,0,0,0) , X(t) = [Bi(t), B 2 {t), Ci(t), C 2 (t), Bi 2 (t), Ci 2 (t)] , (28) 

where vector X describes the part of X that can be obtained by linearization around a classical trajectory. Variations 
near a stable limit cycle rapidly approach zero and therefore X(t) — > 0. As a result, we have only a diffusive growth 
of cumulants Bj (j = 1, 2) as 

Bj(t) =0.5-fjt, (29) 

where we considered the case of a quiet reservoir (rid). This result indicates that the influence of reservoir zero-point 
energy on the dynamics of the nonlinear system is principal physical mechanism responsible for the difference between 
the classical and quantum dynamics in the semiclassical limit. A time scale t* for a correct description of the dynamics 
of quantized SHG in terms of classical electrodynamics can be found by using criterion (|l6|). Taking into account that 
\z(t)\ ~ 1, we have t* ~ 2N-J- 1 . 



Note that the quantum corrections to the classical motion equations (27a) and ( p7b| ) ) do not include cumulants 
B\ i2 . Therefore, in the first order of 1/N, there is no difference between the evolution of quantum mean values and 
the classical dynamics for limit cycle. In other words, the quantum correction Q — > 0, and therefore criterion jl5| ) 
of validity of 1 /TV-expansion does not work. In this respect, the quantized SHG is a somewhat singular problem. In 
other quantum optical systems, for instance, for a nonlinear oscillator with I > 1, typically both criteria of validity 
( |l6| ) and give the same result. 

Over a decade ago, Savage addressed the same problem of quantum-classical correspondence at self-oscillations in 



SHG numerically 14 1. He calculated Q distribution function in the Gaussian approximation centered on a deterministic 
trajectory corresponding to a limit cycle. He worked at a large field and small nonlinearity limits, x/71,2 - * 0, which 
correspond to the classical limit \L^ . It is easy to see that the condition x/71,2 — > is consistent with our condition 
N 3> 1, if one additionally considers the natural condition of a not very strong dissipation in Eqs. (|27|), 71,2/5 ^ 1 
together with g ~ 1 [Eq. (p6|)]. In other words, Savage's small parameter x/l corresponds to our large parameter N 
as xH ~ * N^ 1 / 2 . To establish the difference between the classical and quantum dynamics, the motion equations for 
law-order cumulants were obtained in JTij and solved numerically for particular values of the parameters pi] . Based 
on the results of numerical simulations, Savage concluded that it was a quantum diffusion that was mostly responsible 
for the difference between the classical and quantum dynamics in the semiclassical limit. Moreover, his numerical 
estimate for a characteristic time for the classical description scales as (7/x) 2 > which is in a good agreement with 
our analytical result t* = 2 r y~ 1 N. In summary, our analytical results for the quantum-classical correspondence at 
self-pulsing in SHG are quite consistent with the previous numerical investigation of same problem in |]l4| . 



V. CONCLUSION 



We developed the method of 1 /TV-expansion to consider the nonlinear dynamics and nonclassical properties of light 
in dissipative optical systems in the limit of a large number of photons. The method was applied to the investigation 
of squeezing in higher-order dissipative nonlinear oscillators. We would like to note that our method can also be 
directly applied to an important case of nonclassical states generation in a medium involving competing nonlinearities 

We found a time scale of validity of the l/7V-expansion for a classical description of the dynamics of nonlinear 
optical systems with a simple attractor and a limit cycle. For systems with a simple attractor, this time scale is of 
order unity, and for a limit cycle - proportional to a large N. Qualitatively this result can be understood as follows. 
For time of order unity, the trajectory spirals around a stable stationary point with a small amplitude and therefore in 
virtue of the uncertainty principle, the contribution of quantum corrections to the classical motion equations becomes 
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very important. Unlike the previous case, the oscillations corresponding to a limit cycle are often close to harmonic 
and thus their quantum and classical descriptions can coincide for a fairly long period of time. The basic difference 
between the classical and quantum dynamics in the latter case originates from the influence of reservoir zero-point 
fluctuations, which in our notations are of order of 1/N. This result is in a good agreement with the result of earliest 
numerical simulations of self-oscillations in the quantized second harmonic generation Jl4j . Finally, it should be noted 
that our findings are of a rather general nature and can be applied to the investigations of sclf-oscillations in other 
optical systems, for example, in the optical bistability p3|-p5|. 
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